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Abstract 

The propagation and decay of neutral S-mesons can be described in 
terms of quantum dynamical semigroups; they provide generalized time- 
evolutions that take into account possible non-standard effects leading 
to loss of phase coherence and dissipation. These effects can be fully 
parametrized in terms of six phenomenological constants. A detailed 
analysis of selected S-meson decays shows that present and future dedi- 
cated experiments, both at colliders and S-factories, will be able to put 
stringent bounds on these non-standard parameters. 
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1. INTRODUCTION 



The study of physics related to ^-flavored hadrons appears to be a very promising 
testing ground for various fundamental aspects of the standard model. Indeed, a consider- 
able effort has recently been devoted to the realization of specific experimental programs 
in order to probe with high accuracy the 6-sector of the theory. The so-called 5-factories 
are already collecting data, while experiments at colliders will start running in the near 
future and other will be constructed in the coming years, so that very precise data, in par- 
ticular on CP- violating phenomena, will be available for comparison with the theoretical 
predictions. [1-3] 

The neutral S-meson system is in this respect a imique laboratory: the B^-B^ mixing 
and the induced interference phenomena allow, at least in principle, an accurate analysis 
of very small effects. Besides CP- violating phenomena, these could also include effects in- 
duced by physics beyond the standard model, like those predicted by gran-unified theories, 
supersymmctry and even fundamental dynamics of strings and branes. 

In the following, we shall devote our attention to the study of the non-standard effects 
that are induced in neutral S-meson physics by a suitable generalization of the familiar, 
effective time-evolution of ordinary quantum mechanics. This generalized dynamics takes 
into account possible phenomena leading to irreversibility and dissipation that could affect 
various S-meson observables; we shall discuss in detail to what extent present and future 
dedicated experiments will be able to detect such effects. 

Our approach will be phenomenological in nature; in particular, we shall not need 
to discuss in detail how the generalized time-evolution originates from the fundamental 
dynamics. In fact, the effective description we shall discuss turns out to be largely in- 
dependent from the microscopic phenomena responsible for the appearance of the new 
effects. 

This phenomenological point of view is perfectly consistent: not all dynamics that 
generalize the familiar, unitary time-evolution of quantum mechanics are in fact physically 
acceptable; basic requirements need to be enforced. The new time-evolution should always 
transform, in all physical situations, P-meson states into P-meson states, allow forward in 
time composition law, while increasing the system entropy. As we shall see, these conditions 
uniquely fix the form of the generalized dynamics, that turns out to be parametrized in 
terms of six, new phenomenological constants. In more mathematically precise terms, the 
generalized time-evolution takes the form of a quantum dynamical semigroup. [4-6] 

A physical instance in which evolutions of this type are encountered is given by the 
study of open quantum systems. [4-8] Quite in general, an open system can be considered as 
a subsystem S in interaction with a large environment E. Although the complete system 
S + E evolves in time with the standard unitary operator of quantum mechanics, the 
subdynamics of S alone, obtained by a suitable integration over the environment degrees of 
freedom, is rather involved, showing in general irreversibility and memory effects. However, 
when the interaction between the subsystem and the environment can be considered to be 
weak, the dynamics of S simplifies; it can be represented by linear maps, local in time, 
that take precisely the form of a quantum dynamical semigroup. 

This description is very general and can be applied to model different physical situa- 
tions. It has been originally developed in the framework of quantum optics, [8-10] but it has 
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also been successfully used to study statistical models, [4-6] the interaction of a microsys- 
tem with a measuring apparatus, [11-13] non-standard effects in different interferometric 
experiments involving elementary particles (neutrons, [14] neutrinos [15] and photons [16]).^ 
Furthermore, it provides a general framework for the analysis of dissipative effects in the 
evolution and decay of neutral meson systems. [19, 20, 21-25] (For early investigations, see 
[26].) 

The original physical motivation for such investigations came from quantum gravity: 
due to the quantum fluctuation of the gravitational field and the presence of virtual black 
holes, space-time should lose its continuum character at Planck's scale, and this could 
lead to loss of quantum coherence. [26-33] From a more fundamental point of view, also 
the dynamics of strings and branes could act as an effective environment, inducing non- 
standard, dissipative effects at low energies. [34, 35] Fortunately, as already mentioned, the 
details of the microscopic, "stringy" dynamics are not needed for the effective description of 
the dissipative phenomena in terms of quantum dynamical semigroups. [35] In this respect, 
this phenomenological approach provides a universal framework for the study of quantum 
decoherence effects. 

In the case of neutral kaon system, the analysis of the dissipative phenomena has been 
pursued in detail, and upper bounds on some of the constants parametrizing the new effects 
have been obtained using available experimental data. [21, 22] Improvements of these results 
are expected from the experimental study of correlated K-mesons at (/)-factories.[23, 24] 

For the case of the neutral S-meson system, a preliminary analysis of the relevance of 
dissipative effects on semileptonic decays has been presented in [25]. There, a particular 
phase choice for the B states has been adopted; it allowed treating indirect CP- violating 
effects on the same footing as the decoherence phenomena. 

In the following, a much more general and complete discussion will be presented. Using 
a manifestly phase-invariant formalism that makes no commitment on the magnitude of 
CP (and CPT) violating effects, the impact of the irreversible, dissipative phenomena on 
neutral S-meson decays will be analyzed and discussed in detail, using both semileptonic 
and selected hadronic decay channels. We shall first study decays of single mesons, relevant 
for experiments at colliders, and then discuss the case of correlated mesons at 5-factories. 
In both cases, the new, dissipative phenomena modify in a specific and characteristic 
way the various S-meson observables, so that the presence of dissipation can be probed 
quite independently from other non-standard effects. -l- Although more specific studies, 
that include precise analysis of acceptance and efficiency of the various detectors, are 
certainly necessary, our investigation clearly indicates that present and future S-meson 
experiments should be able to put stringent limits on the parameters that describe non- 
standard, dissipative phenomena. 



T Propagation of neutrinos in thermal environments can also be described by equations 
of semigroup form; for details, see [17, 18] and the discussion in [15] 

* For recent studies on possible violations of the CPT symmetry and other, unconven- 
tional phenomena in S-decays, see [36-41] and references therein. 
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2. QUANTUM DYNAMICAL SEMIGROUPS 



The familiar effective description of the propagation and decay of neutral 5-mesons 
requires the introduction of a two-dimensional Hilbert space; ^ the time-evolution is then 
realized by a linear transformation on the elements of this space. [1, 2] In presence of 
dissipation however, a more general formalism is needed, in which S-meson states are 
represented by density matrices. These are hermitian, positive operators {i.e. with positive 
eigenvalues), with constant trace (at least for unitary evolutions). In the \B^) basis, 

any S-meson state can then be written as 



Pi P3 

Pa P2 



(2.1) 



where p4 = pg, and * signifies complex conjugation. 

As explained in the introductory remarks, our analysis is based on the hypothesis that 
the time evolution of the state p is realized in terms of linear transformations: they are 
generated by an equation that extends the familiar quantum mechanical one, being of the 
form: 

dp{t) 



dt 



-iH p{t) + ip{t) + L[p] . (2.2) 



The first two pieces on the r.h.s. give the standard hamiltonian contribution, while L is a 
linear map that encodes possible dissipative, non-standard effects. 
The effective hamiltonian H includes a non-hermitian part, 

H = M-'-r , (2.3) 

with M and F hermitian matrices, that characterizes the natural width of the states. 
The entries of H can be expressed in terms of its eigenvalues: Xs = rns — fTs, Xl — 
niL — |7l, and the complex parameters ps, Qs, Pl, Ql, appearing in the corresponding 
(right) eigenstates, 

\Bs) = PS 1^0) +qs\B^) , |P5|' + = 1 , 

|2 ^ ^ ^^-^^ 



\BL)=PL\B'')-qL\B^) , bLl' + k^l 



We use here a notation that follows the conventions usually adopted for the neutral kaon 
system. The two states in (2.4) are expected to have a negligible width difference, 

Ar«r, Ar = 75-7L, y = ^^^^, (2.5) 

so that they may be more conveniently distinguished via their mass difference. Am = 
rriL — ms, rather than their different lifetimes. For this reason, the notation for the 

't' We shall limit our considerations to S^- mesons, although most of the discussions 
below could be applied to Sg-mesons as well. 
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light meson, and \Bh)i for the heavy partner are therefore sometimes used instead of those 
in (2.4). [1, 2] 

The effective hamiltonian H can be diagonahzed using the similarity transformation 
induced by (2.4): 

(2.6) 
(2.7) 



H = VHqV 



-1 



with 

Then, one can write 
H = 



V 



PS PL 

qs -ql 



Hn 









Xl 







1 


rsXs + tlXl 


rsr^Xs 


-Xl) 




H2 


rs+rL 


As — Xl 


tlXs + 


tsXl 



where we find useful to introduce the two ratios: 



PS 

rs ^ — 
Qs 



PL 

rL = — 
Ql 



(2.8) 



(2.9) 



They can be conveniently expressed as 



rs 



1 + e 

1-9 



tl = vo- 



i-e 



(2.10) 



in terms of the two complex parameters 



rsTL , 



9 



rs - tl 



rs + rL 

that signal T and CPT violating effects. Indeed, from (2.8) one has: 



(2.11) 



^3 f> 



Hi — H2 
Xs — Xl ' 



(2.12) 



so that T (and CP) invariance is broken when \a\ ^ 1, while CPT (and CP) invariance 
is lost for ^ 7^ 0.[1, 2] 

Even in absence of the additional piece L[p] in (2.2), probability is not conserved during 
the time evolution: dTT[p{t)]/dt < 0. This is due to the presence of a non-hermitian part 
in the effective hamiltonian H. On the other hand, loss of phase coherence shows up only 
when the piece L[p] is nonvanishing: it produces dissipation and possible transitions from 
pure states to mixed states. 

As already mentioned before, not all maps L[p] gives rise to integrated time evolutions, 
7t : p(0) 1-^ p{t), that are physically acceptable. Quite in general, the one parameter 
(=time) family of linear maps 7^ should transform S-meson states into S-meson states, 
and therefore should map an initial density matrix into a density matrix; further, it should 
have the property of increasing the (von Neumann) entropy, S = — Tr[p(t) Inp(t)], of 
obeying the semigroup composition law, 7t[p(t')] = p{t + t')i for t, t' > 0, of preserving the 



5 



positivity of p{t) for all times. Actually, for the physical consistency of the formalism in 
the case of correlated systems, one has to demand that the time evolution 7^ be completely 
positive. [4-6] Once these properties are taken into account, the form of L[p\ results uniquely 
fixed and the family of maps 7^ represents a so-called quantum dynamical semigroup. 

The linear map L[p] can be fully parametrized in terms of six real phenomenolog- 
ical constants, a, b, c, a, j3 and 7, with a, a, 7 non negative, satisfying the following 
inequalities: [21, 23] 

2i? = a + 7- a>0, RS-ll^yQ , 

25 = a + 7- a>0, RT - > Q , 

" (2.13) 
2T = a + a-7>0, ST - {3^ >Q , ^ ' 

EST -2hcj3- R(3^ - Sc^ -Tb^>0 , 



direct consequence of the property of complete positivity. A convenient explicit expression 
for L[p] can be obtained by introducing a vector notation for the matrix p, rewritten as a 
four-dimensional column vector \p), with components (pi, P2, Ps, Pa)- In this way, the map 
L[p\ can be represented by a 4 x 4 hermitian matrix C, acting on \p): 



—a 
a 

-{c + ib) 
l-{c-ib) 



a 
—a 
c + ib 
c — ib 



ib) 
c — ib 



c 

-(a + 7) 
a-^ + 2ip 



a 



ic + ib) ■ 

c + ib 
- 7 - 2i/3 

-(a + 7) ■ 



(2.14) 



Among the physical requirements that the complete time evolution 7^ should satisfy, 
complete positivity is perhaps the less intuitive. It is often dismissed, in favor of the more 
obvious simple positivity. Simple positivity is in fact generally enough to guarantee that 
the eigenvalues of the density matrix p{t) remain positive at any time; this requirement is 
obviously crucial for the consistency of the formalism, in view of the interpretation of the 
eigenvalues of p{t) as probabilities. 

Complete positivity is a stronger property, in the sense that it assures the positivity of 
the density matrix describing the states of a larger system, involving the coupling with an 
extra, auxiliary finite-dimensional system. Although trivially satisfied by standard quan- 
tum mechanical (unitary) time-evolutions, the requirement of complete positivity seems at 
first a mere technical complication. Nevertheless, it turns out to be essential in properly 
treating correlated systems, like the two S-meson coming from the T(45') resonance; [42] 
it assures the absence of unphysical effiects, like the appearance of negative probabilities, 
that could occur for just simply positive dynamics (see the discussion in Sect. 7). 

A further comment on the expression (2.14) of the dissipative part of the evolution 
equation (2.2) is in order. The entries of the matrix (2.14) have been written in a specific 
choice of basis in the two-dimensional Hilbert space, the one for which the state p takes the 
form (2.1). If one performs a (unitary) transformation on the basis vectors \B^), \B'^), also 
the six parameters appearing in (2.14) will in general change to new ones a', b', c', a', /?', 
and 7', expressed as linear combinations of the old ones. Despite this, one can check that 
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the inequalities in (2.13) are form- invariant; in other terms, the transformed dissipative 
parameters still obey (2.13), provided the old ones do.^ 

This discussion clearly illustrates that, although the equation (2.2) has been written in 
a fixed, specific S-meson basis, the notion of complete positivity is basis-independent, as is 
the description of the dissipative effects: if the contribution L in (2.14) is non vanishing in 
one basis, it is non vanishing in any basis. This result has to be expected from the general 
theory of quantum dynamical semigroups, that can be formulated in a basis- independent 
way, [4-6] and will be made apparent in the discussion of the coming sections. 



3. THE EFFECTIVE TIME-EVOLUTION 

The behaviour in time of physical observables related to the various S-meson decay 
channels can be obtained by solving the evolution equation (2.2) for an arbitrary initial 
state p(0). This results in the study of a system of linear differential equations for the 
entries of the density matrix p in (2.1). 

It is convenient to use the vector notation introduced in the previous section, and 
write the matrix p as the four-dimensional vector |p). Then, the evolution equation (2.2) 
takes the form of a Schrodinger (or diffusion) equation: 

||p(t))=/C|p(t))^[7i + £]|p(t)), (3.1) 



where H. is the 4x4 matrix containing the hamiltonian contributions [c/. (2.8)], 

(3.2) 



n = 



2Im{Hi) iH^ -iHs 

2Xm{H2) -iH^ iHl 

iHl -iHs i{H^-Hi) 

L -iH4 iHl i{H*-H2)] 



while the dissipative part C is given in (2.14). 

As already mentioned before, the hamiltonian piece in (3.2) contains contributions 
that are not invariant under CPT and T transformations; this is also true for £, so that 
in general dissipation will induce violations of these discrete symmetries. It is instructive 
to explicitly discuss this point in the formalism of (3.1).* 

Following the rules of quantum mechanics, any symmetry transformation can be re- 
alized by a unitary or antiunitary operator U acting on the basis states \B^), l-B*^), or 
alternatively on the density matrix p. This induces an action on the vector |p), realized 
by a 4 X 4 matrix lA: 

|p) ^ |p') =U\p). (3.3) 



t For a generic, non-unitary change of basis, as the one discussed in Appendix A, also 
the form of the conditions (2.13) in general changes, while always assuring the fulfillment 
of the property of complete positivity. 

^ Further discussions on the notion of symmetry invariance for quantum dynamical 
semigroups can be found in [43]. 
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This transformation leaves the equation (3.1) form-invariant, i.e. it is a symmetry of the 
evolution equation, provided: 



for unitary transformations, or alternatively 

in the case of antiunitary transformations. 

For instance, an independent phase change of the two basis vectors. 



(3.4a) 
(3.46) 



is realized by the 4x4 diagonal matrix: 



|S0) ^ e^-^lSO) , 



U^, = diag 



1, 1, e^^'^~'^\ Q-^i4>-<t>) 



(3.5) 



(3.6) 



One easily checks that this transformation is not an invariance of Ti and hence of (3.1); 
indeed, the off-diagonal elements and of the effective hamiltonian H do change 
under (3.5). As a consequence neither cr in (2.12), nor rs and tl are phase invariant. 

In the case of discrete CPT, T and CP transformations, one finds that the corre- 
sponding matrix U takes a block diagonal form: explicitly, one finds: 



UcPT = 



fTi 
(To 



Ut = 



(Jo 



Ucp = 



(71 



(3.7) 



where cXi, i — 1, 2, 3 represent the usual Pauli matrices, (Tq being the identity, while 



(Tip 



11 e-^'^f 

f.2i^ 



(3.8) 



contains the dependence on the phase that defines the CP transformation of the B^-B^ 
basis states: ^ e^^^^B^). 

For the hamiltonian contribution 7i, insertion of the expressions (3.7) in the appro- 
priate invariance condition (3.4) gives the familiar results. In particular, CPT invariance 
requires Hi = H2, or equivalently, recalling the definitions (2.11), (2.12): 



rs = rL , 

while T invariance implies {H^l = {H^l, i.e. 



9=0 



\rsrL\ = 1 , 



I 1-1 + ^ 



1 







(3.9) 



(3.10) 



Further, notice that imposing both CPT and T invariance, hence CP invariance, i.e. 
= ^ = 0, readily implies: 



rs = rL = , 



(3.11) 
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For the dissipative part C, of the pseudo hamiltonian /C similar invariant conditions can 
be derived using (3.4a) (the matrix in (2.14) is in fact hermitian); they impose constraints 
on the six dissipative parameters a, 6, c, a, (5 and 7. More specifically, invariance under 
CPT-transformations requires the vanishing of c and 6, while T-invariance imposes the 
conditions: (c - ih) = e^''^{c + ih) and (a - 7 + 2i(3) = e^^'^(a - 7 - 2i/3). Finally, CP- 
invariance requires (c — i6) = — 6^**^(0 + z6) and (a — 7 + 2i/3) = e'^*'^(Q; — 7 — 2i/?). It should 
be noticed that all these constraints are perfectly compatible with the inequalities in (2.13); 
in other words, complete positivity does not interfere with the discrete transformations. 

Although the available experimental bounds on the CPT and T violating contributions 
of the hamiltonian H are not very accurate, the magnitude of the constants 9 and ^ 
parametrizing these violations are expected to be very small. Despite this, in the discussion 
that follows we shall try to be as general as possible, and keep 9 and ^, hence rg and r^, 
arbitrary, unless explicitly stated. 

In solving the evolution equation (3.1) however, we shall assume the dissipative con- 
tribution C to be small. In a phenomenological approach, it is hard to give an apriori 
estimate on how large the dissipative effects should be. However, as already mentioned 
in the introductory remarks, a general framework in which dissipation naturally emerges 
is provided by the study of subsystems in interaction with large environments. In such 
instances, the non-standard effects can be roughly estimated to be proportional to powers 
of the typical energy of the system, while suppressed by inverse powers of the characteristic 
energy scale of the environment. [4-7] 

In the case of the B^-B^ system, these considerations, together with the general 
idea that dissipation is induced by quantum effects at a large, fundamental scale Mp, 
lead to predict very small values for the parameters a, b, c, a, (5 and 7; using dimensional 
arguments, an upper bound on the magnitude of these parameters can be roughly evaluated 
to be of order m\ /Mp^ where ms is the neutral S-meson mass. If the non-standard effects 
have a gravitational origin, the scale Mp should coincide with the Planck mass Mp, and 
the previous upper bound would give: m\/Mp ~ 10~^^ GeV. 

These considerations allow treating the dissipative piece C in (3.1) as a perturbation to 
the hamiltonian contribution Ti. In order to set up the perturbative expansion, it is useful 
to make a change of basis in (3.1), so that the hamiltonian piece Ti becomes diagonal. 

As discussed in the previous section, the effective hamiltonian H can be brought to 
diagonal form by the similarity transformation (2.6); the diagonalizing matrix V in (2.7) 
can be conveniently decomposed as: 



V = V-Q, V = 



rs rp 
1 -1 



Q = 



qs 
QL 



(3.12) 



Since the matrix Q is diagonal, it disappears from the relation (2.6), that can be equally 
well be written as: 

H^VHoV-^. (3.13) 

When applied to the density matrix p, this change of basis induces the transformation 
p ^ p — p or equivalently: 

|p)^|p) = V|p) , (3.14) 
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where 



V 





-1 


1 1 2 

\rL\ 


rl 




1 


1 


1 |2 


-r*s 


-rs 


\rs + rL\'^ 


1 


-r%rL 




tl 




.1 


-rsrl 




-rs- 



(3.15) 



Then, the evolution equation (3.1) becomes 

d 



dr 



\~p{t))= no + c\\~p{t)) 



with 



(3.16) 
(3.17) 



and a diagonal hamiltonian piece: 



-is 








- 





-IL 














-r_ 














-r+J 



r±zAm 



Am = rriL — ms 



(3.18) 



(3.19) 



The price to pay for this change of basis is a more complicated expression for the matrix 
jC, representing the dissipative contribution (its explicit form is collected in Appendix A). 
However, the entries of C appear now to be manifestly invariant under the change of phase 
in (3.5), and therefore, the same is true for the whole evolution equation in (3.16). This 
is a great advantage since any solution of (3.16), even an approximate one obtained using 
perturbative methods, will result to be manifestly phase-invariant. 

We remark that physical observables, being the result of a trace operation (see below), 
are by definition independent from any phase convention. Provided a sufficient number 
of observables are computed, one can consistently re-express a result obtained in a given 
phase convention, into any other phase choice. However, once a phase choice is adopted 
and a given approximation used, conclusions drawn from only a limited number of physical 
observations could lead to incorrect conclusions. It is therefore always preferable at any 
step to work in a phase-independent framework. 

The evolution equation in (3.16) can now be solved by iteration, to any order in the 
small dissipative parameters a, 6, c, a, (3 and 7. Simple manipulations allow to express 
the vector \p{t)) at time t as the series expansion: 



|p(t))=e^o* |p(0))+ / dse^°^'-''^ JZe^^'' |p(0)) 

Jo 

+ [ dsi r ds2e^°^'-''^ jCe^°^''-''^ jCe^°'' |p(0)) + ... . 
Jo Jo 



(3.20) 



For the considerations that follow, it will be sufficient to consider terms that are at most 
linear in the small parameters, thus retaining only the first two terms in the r.h.s. of (3.20). 
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This suggests a further simphfying assumption, that allows to write down a more manage- 
able expression for the matrix L. 

As it is apparent from its definition in (3.17), the entries of L are linear combinations 
of a, 6, c, Q!, /3 and 7, with coefficients that depend on 6*, ^ and the phase of a . Since Q 
and ^ are themselves expected to be small, they can be safely ignored in those expressions. 
In other terms, it appears reasonable to neglect CPT and T (hence CP) violating effects 
induced by non- vanishing d and ^, when these violating phenomena mix with the dissipative 
effects. 

Within this approximation, the matrix L takes the form 



-D 


D 


-c 


-C* 


D 


-D 


c 


C* 


-C* 


C* 


-A 


B 


-c 


c 


B* 


-A 



(3.21) 



where 

1 = i I [(a + 7) + 2a + 7^e[(a - 7 + 2i^)a\ | , (3.22a) 

5 = ^ I [(a + 7) - 2a + ■Re[{a - 7 + 2i/3)(j] - 4iXm[(c - %h)^Ja\^ , (3.226) 

C ^ ^|2 7le[(c-i6)ya] + iXm[(a-7 + 2i/3)a]| , (3.22c) 

5 = ^ I (a + 7) - -ReXipi - 7 + } , (3.22d) 

and cr is here the pure phase defined in (3.11). 

It is now a matter of a simple computation to perform the integrals in (3.20) and 
find the time dependence of the components P\{f), ^2(^)7 ^3(^)7 P4(t) of the vector |p(t)), 
that give the entries of the density matrix p(t). The explicit expressions are collected 
in Appendix B, and will be used in the next sections to study in detail various neutral 
S-meson decays. 



4. OBSERVABLES 

In the formalism of density matrices, any physical observable of the neutral S-meson 
system is described by a suitable hermitian operator O. Its evolution in time can be 
obtained by taking the trace with the density matrix p(t). 

Of particular interest are those observables Of that are associated with the decay of 
a S-meson into final states /. In the \BS')^ basis, O/ is represented by a 2 x 2 matrix. 



O4 O2 



(4.1) 



11 



whose entries can be explicitly written in terms of the two independent decay amplitudes 
A{B^^ f) and A(B^ ^ f): 



02 = \A{W ^ , A{B' - /) [A{W ^ /)] * . 

Being a physical quantity, directly accessible to the experiment, its mean value, 

{Of)^Tr[Ofp] , 



(4.2) 



(4.3) 



is however basis independent; it can be computed in any specific representation. In par- 
ticular, using the transformation V introduced in (3.12) in the definition (4.3), the time 
evolution of (Of) can be written as 



4 



(4.4) 



i=l 



where the entries Oi of the transformed matrix Of — V^OfV have been labelled as in 
(4.1). This general formula will be repeatedly used in the coming sections to explicitly 
compute experimentally relevant decay rates and asymmetries. 

In parametrizing B^-B^ decays, it is customary to introduce the following phase- 
independent complex quantities: [1, 2] 



f ^ qs AjBO /) 
PS ^(50 ^ /) ' 



A 



J _ A{BO f) 



^ PL A{B^ ^ /) • 
They can be used to express the entries of the matrix C/; explicitly, one finds: 



(4.5) 



Of = \A{B'^f)\ 



\fs\ 



|i + a{| 



r%rL (1 + A^)*(1-A{) 



rsrl (1 + A^) (1-A{)* 
It is sometimes convenient to use instead the parameters: 



rL\ 



(4.6) 



f ^ PS A{B^ - 
qs A(B^ - 

one can then rewrite (4.6) as: 

df = \A{W^f)\' 



f) 



f) 



f PL AiB^' 

f^L= 



/) 



1 1 f |2 

|1+4| 



QL A{B^ ^ /) Xi 
h / 12 



(4.7) 



(4.8) 



The form (4.6) and (4.8) for the observable Of are very general and hold for any final 
decay state /. Simplified expressions can however be derived by looking at specific decay 
channels. 
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Let us first consider decays of the neutral S-mesons in semileptonic states Mv^ where 
h stands for any allowed charged hadronic state. We shall be as general as possible, and 
include in our discussion violations of the AS = AQ rule. The amplitudes for the decay 
of a or B'^ state into h~£^iy and h^£~i' can then be parametrized in terms of three 
complex constants, Xh, Vh ctnd Zh, as follows: 

Mh{l-yh), (4.9a) 

MUl^yl) , (4.96) 

ZhA{BO ^ h+e-u) , (4.9c) 

XhA{B^ ^ h-£+u) , (4.9d) 

where A4h is a common factor. (Note that sometimes, [44, 1, 2] the notation Xh = zj^ 
is used instead of Zh-) The AS = AQ rule would forbid the decays B^ h'^i~u and 
B'^ h~£~^u, so that the parameters Xh and Zh measure the violations of this rule. 
Instead, CPT-invariance in the decay process would require i/h — 0. In view of this, the 
quantities Xh, Uh ^md Zh are expected to be very small. In the computation of semileptonic 
observables, it seems therefore justified to keep only first order terms in these parameters, 
neglecting terms that contain Xh, Vh and zu multiplied by the dissipative parameters 
or the constants 9 and ^, that signal CPT and T violations in "mixing". Within this 
approximation, the quantities in (4.5) and (4.7) reduce to 

><f = Al" = V^Xh = Xh , l^t = l^f = V^^h = l^h , (4.10) 

where a is the pure phase of (3.11). 

For hadronic final states, such simplifications are in general not possible. Neverthe- 
less, to be consistent with the approximation adopted in the previous section, when the 
coefficients A5, Al, //g and hl multiply a dissipative parameter, they should be computed 
in the limit of exact CPT and CP symmetries. This is particularly relevant when the final 
state / has a definite CP-parity Q; within that approximation, one in fact finds: [1] 

Xi = Xl = Cf. (4.11) 

As we will see in the following, this observation will be helpful in the explicit evaluation 
of the asymmetries involving decays into these specific final states. 



A{B^ h-£+u) = 
A{W ^ h+£-u) = 
A{B^ h+£-u) = 
AlW h-£+u) = 



5. SINGLE MESON DECAYS 

In this section, we shall study observables connected to the time-evolution and decay 
of a single, uncorrelated B^-B'^ system, that can be typically measured at colliders. Let us 
indicate with \pBo{t)), \pBo{t)) the time evolution according to (3.16) of the states Ip^o), 
Ip^o) that represent initial pure B^, B^ mesons; their components can be organized in the 
two matrices: 

1 

Pbo = I ^ 

\rs + tl\ 



1 1 
1 1 



\rs + rA 



-rsrl 



(5.1) 
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As discussed in the previous section, the probabihty rate that an initial \pb°)j \Pb°) state 
decays at time t into a given final state / described by the operator Of is given by: 



Vf{B';t)=Tr[OfPBo{t)] , VfiBO;t)=TT[Ofpso{t)] 



(5.2) 



In writing down the explicit form of the probabilities V, it is convenient to introduce the 
new variable r = tF; in practice, r expresses the time variable in units of the B lifetime. 
It is also customary to define the two combinations:''' 



Am 



u> 



AT 
2F 



(5.3) 



Although not yet directly measured, d is expected to be very small, d < 10~^. For this 
reason, in discussing B^-B^ observables, one often takes the simplified assumption d — 0. 
This choice produces however important consequences. Indeed, in the limit AF = 0, 
the antihermitian part F of the effective hamiltonian H in (2.3) becomes proportional to 
the identity. In this case, the similarity transformation (2.6) that diagonalizes H turns 
out to be unitary and the eigenstates \Bs) and \Bl) in (2.4) become orthogonal. Since, 
{Bl\Bs) = Ps Pl '~ Qs Qlj ^^^^ implies the condition 



* -1 



(5.4) 



and as a consequence, |cr| = 1, or equivalently ^ = 0, and also Xm{6) = 0. Therefore, 
the condition 5=0 implies T conservation in the hamiltonian piece of the evolution 
equation, while CPT violation is possible only if TZe{9) is different from zero. This fact 
has clearly important consequences for experimental tests of CPT and T invariance in 
"mixing" .[45, 46] In the following, we shall keep 6 nonvanishing, unless explicitly stated. 

We shall first consider observables connected with S-meson decays into semileptonic 
final states. Using the appropriate observables (4.6), (4.8) and the results in (4.10) together 
with those collected in Appendix B, expressions for the probabilities (5.2) can be explic- 
itly obtained. Apart from a common exponential decay factor, they show an oscillatory 
behaviour modulated by co and further exponential terms regulated by S: 



cos UJT 

+ sin u)T 
+ cosh 5t 
+ sinh 5t 



\fs -\-tl\ I , .2 



52 + 0;^ 



4Tm(r«n) _^ 



2(|rs 


2 + 








rs + tl 


2 



2ne{yh) + 



46 



,ne{C) 



2(|r5 


2 _ 


rL 






rs +rL 


2 



-2ne{Xh) + 



D 



(5.5a) 



't' They are usually called xb and j/s, respectively; we do not use these labels to avoid 
confusion with the parameters introduced in (4.9). 



14 



T^h+ [B ; r) = — - — e 



cos LOT 

+ sin COT 
+ cosh 5r 
+ sinh 5t 



^^-{A-D)t 

Ws + rA'^ 



2Im{^ih) - ne{B) + 



4.5 



52 + a;2 



Xm(C) 



+ 2ne{yh) 



0^ + a;^ 



(5.56) 



where the effects of the non-standard, dissipative phenomena are controlled by the dimen- 
sionless parameters [cf. (3.22)] 



B 



B 

Am 



C 



C 



D 



D 



(5.6) 



The expressions for Vh-{B^;T) and Vh+{B^;T) are obtained from (5.56) and (5.5a), re- 
spectively, by changing the sign of yh and C, and letting rs — > 1/^5, tl — ^/tl and 
A/i ^ Hh- One can check that the formulas (5.5) reduce to those presented in [25] when 
adopting the particular phase convention for which 



rs ~ 1 + 2 €5 , ri, ~ 1 + 2 , 



(5.7) 



and the quantities es and ez,, parametrizing CPT and CP violation in "mixing", are 
considered small. Instead, we stress that the expressions in (5.5) are completely general, 
and manifestly independent of any phase-choice. 

The probabilities Vh are directly accessible and can be studied in experiments per- 
formed at colliders. A preliminary investigation, assuming no dissipation, has already been 
performed at LEP in the framework of the approximation (5.7) and an upper bound on the 
magnitude of 7le{es + cl) has been obtained. [47, 48] Much more accurate studies can be 
performed both at the Tevatron and LHC, where the r-dependence in (5.5) can actually 
be observed; extrapolating from preliminary simulation estimates, sensitivities of a few 
percent on at least some of the parameters in (5.6) can be reasonably expected. [49, 50] 

Instead of analyzing directly the decay rates in (5.5), it might be more convenient 
to study specific asymmetries, constructed by taking suitable combinations of the VhS. 
Two independent simple asymmetries that can be formed with the probabilities Vh and 
are likely to be studied at colliders are: 



^cpt(t) = 



P^+(SO;t)+P^-(SO;t) 



Vh-{B';r)+VhAB';r) 



(5.8a) 



(5.86) 
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The first one is sensible to CPT violating effects, by comparing the rate difference between 
the process and its CPT-conjugate B'^ B^, while the second asymmetry 

signals the violation of time reversal by measuring the rate difference between the process 
B^ — > B^ and its time-conjugate B^ — > B^. 

Unfortunately, these simple observables are not very sensitive to the presence of the 
non-standard, dissipative effects. For instance, by using the expressions in (5.5), one finds: 





2 _ 


rL 


2) sinh 5r + 2 Xm(r5 r2) sinwr 


{\rs 


2 + 




2) cosher + 2 7^e(r5r2) coswr 



acpt{t) = , , + ^neiyu) 

L 

1 



+ 



cosh 5t + cos LOT 

( cosh (5t — COS a;T \ 45 



\cosh 8t + cos LOT ) + up' 



Xm{Xh — f^h) sina;T + TZe{Xh — sinh^r 
ne{C) , 



(5.9) 



which depends on the dissipative parameters through TZe{C). However, in (5.9) this pa- 
rameter appears multiplied by a factor 5, so that its contribution to Acpt{j) turns out 
to be suppressed with respect to the other terms; in particular, the approximation 5 ~ 
would completely eliminate the presence of this parameter from (5.9). The same conclusion 
holds for the other asymmetry At- Therefore, observables of type (5.8) are not suitable 
for probing the presence of dissipative effects in neutral 5-meson physics. 

More complicated combinations of the probability rates Vh are needed in order to 
isolate the dissipative contributions. A particularly illuminating example is given by: 

^ . . ^ [Vh- (i?°; r) - V,- (BO^r)] - [V,, {B^; r) - P,. (g^; r)] 
^"^^^^ P^-(SO;r) + n-(SO;r)+P,+ (SO;T) + n+(SO;T) 

Even in the approximation S = 0, this asymmetry retains a distinctive dependence on the 
parameters that signal the presence of dissipation: 

AAmir) ^ [ne{e)Y {1 -cosiOT)+e ^'^ cos LOT + [JZe{B) - I'm{\h + iJ'h)] sin uur . (5.11) 

Assuming the validity of the AB = AQ rule, and neglecting the quadratic dependence 
on the CPT- violating parameter 9, the different time behaviours in AAm{T) should allow 
a determination of the non-standard parameters A and the real part of B. The actual 
accuracy of such a determination highly depends on the sensitivity of the measure of ^Am- 
Although specific studies have not yet been performed for the form (5.11) of this asym- 
metry, from the results of available simulations one can nevertheless expect an accuracy 
in the determination of the various terms in Aa™, of about a few percent. [49, 50] This 
sensitivity is already enough to give interesting bounds on A and 7le{B). 

Nevertheless, the most simple and accurate tests on the extended dynamics in (3.1) 
that can be performed at colliders, using semileptonic decays, is based on the analysis of 
time-integrated rates, defined in general as 

'PAB)^^ dTVf{B-T) . (5.12) 
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Despite most asymmetries constructed with these integrated probabihties suffer the same 
(5-suppression problem mentioned before in discussing the observables in (5.8), a combi- 
nation analogous to the one in (5.10) is again very sensible to the presence of dissipative 
parameters. An explicit evaluation gives, for 6—0: 



A\ - 

■^Am — 



+ 



KeiB) - Xm{\h + iJih) 



{uj^ -l)A-2uj^D 



(5.13) 



1 + 

Assuming again the AB = AQ rule and neglecting quadratic CPT- violating effects in 
"mixing" {[JZe{6)]'^ = 0), a measure of A'^^ (1 + a;^) not compatible with 1 would clearly 
signal the presence of non-standard effects. Such a test can be easily performed at colliders 
with high accuracy, in particular using the next generation dedicated S-meson experiments. 

Further useful information on the dissipative parameters in (5.6) can be obtained 
by studying observables and asymmetries involving S-meson decays into hadronic final 
states /. The explicit form of such observables are in general more involved than the ones 
containing semileptonic decays, since now one can not resort to simplifying approximations 
as in (4.10). Nevertheless, when the final states / exhibit a definite CP-parity Q, at least 
within a certain approximation, manageable expressions for the relevant asymmetries can 
be obtained. 

The typical observable that can be constructed in this case involves the difference of 
the two probabilities in (5.2): 



Mr) = 



(5.14) 



Vf{BO;T) + Vf{BO;T) 

where, for instance, / can represent D'^D~, tt+tt" or J/ifjK final states. These asym- 
metries will be the target of very intensive experimental studies: in fact, they offer the 
possibility of clean tests of the standard model paradigm for CP-violation. On top of this, 
they turn out to be very sensible to the determination of the dissipative parameter C in 
(5.6). 

For simplicity, in presenting the explicit expression for the asymmetry A/(r) we shall 
take S = and retain only first order terms in the CPT-violating parameter TZe{9). Within 
this approximation the two quantities A;^ and A{ parametrizing the B^-B'-^ decay into / 
can be expressed in terms of a unique amplitude ratio, 

- A{W f) 



as follows 



Xf [1 - TZeiO)] , Xi = Xf[l + TZeid)] 



(5.15) 



(5.16) 



Using these results in the computation of the probabilities in (5.2), one finally obtains: 



A fir) 



2^e(A^) 
1+|A/|^ 



ne{e)+^Im{C) - 



2Im{X f) , 2Cf 



+ 



a; 



ne{C) 



sm UT 



+ 



1-lA/p 2ne{Xf) 
1 + lA/P l+\Xf\^ 



TZe{e) 



2C/ 



(jj 



Xm{C) 



(5.17) 



cos ijJT 
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The presence of dissipation manifests itself in the modification of the coefficients of the 
oscillating terms and in the presence of a r-independent piece. Neglecting CPT- violations 
in "mixing" {6 = 0) and assmning a vanishingly small "direct" CP-violation, i.e. |A/| ~ 1, 
the different r-dependence in the three terms of A fir) allows a determination of both 
TZe{C) and Xm{C), together with the CP-violating parameter Xm{\f). Dedicated B- 
experiments at colliders should be able to identify the various r-dependence and therefore 
provide very stringent bounds on the parameter C. 

On the other hand, integrated asymmetries, constructed from the probabilities in 
(5.12), are not very useful in probing the presence of dissipative effects in neutral P-meson 
decays. For instance, using the same approximations that lead to (5.17), the integrated 
asymmetry A'^ defined as in (5.14) takes the form: 



1 - 


A/ 


2 


1 + 


A/ 


2 



Even for vanishing TZe{9) and |A/| = 1, the measure of this asymmetry alone gives little 
information on the magnitude of C. 



6. CORRELATED MESON DECAYS 



The non-standard, dissipative effects described by the generalized dynamics (3.1) can 
be further analyzed in experiments involving correlated B^-B^ mesons, at the so-called 
P-factories. Indeed, as we shall see, these set-ups appear particularly suitable for studying 
phenomena involving loss of quantum coherence. 

In those experiments, correlated P°-PO mesons are produced from the decay of the 
T(45') resonance. Since the T-meson has spin 1, its decay into two spinless bosons pro- 
duces an antisymmetric spatial state. In the T-rest frame, the two neutral P-mesons are 
produced flying apart with opposite momenta; in the basis |P°), |P'^), the resulting state 
can be described by: 



\iPa) = ^ (l^o, -p) ® |po,p) - |P0, -p) ® |pO,p)) 



(6.1) 



The corresponding density operator pA can be abstractly expressed in terms of a projector: 



PA = IV'a) (V'aI 



(6.2) 



As in the previous sections, we find it more practical to use the transformation introduced 
in (3.12), (3.13) and pass to a different representation: 



PA 



PA = 



V 



PA 



(6.3) 
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The 4x4 matrix pA can then be exphcitly expressed as: 

1 



PA = 



2 \rs +rLp 



Pi® Pa + P2®Pi - Ps^Pa - Pa®Pz 



where 



Pi 



Ps = 



1 



1 





P4 = 





1 



1 



(6.4) 



(6.5a) 



(6.56) 



The evolution in time of the density matrix can be analyzed using the single B- 
meson dynamics p{t) discussed in Sect. 3. We shall assume that once produced in a T-decay, 
the two mesons evolve in time each according to the completely positive map 7t generated 
by (3.16).^ This guarantees that the resulting dynamics is completely positive and of 
semigroup type. As already mentioned, this assures the positivity of the eigenvalues of 
any physical states at all times: this is an essential condition for the consistent description 
of any physical system. [23, 42] 

The density matrix that describes a situation in which the first P-meson has evolved 
up to proper time ti and the second up to proper time t2 is then given by: 



pAitiM) = (7ti ® Tta) [pa] 
1 



2|rs + rL|2 



Pl{tl) ® P2{t2) + ^^2(^1) ® Pl{t2) 

- Psih) ® P4(t2) - hih) ® P3(t2: 



(6.6) 



where Pi(ti) and Pi{t2), i = 1,2,3,4, represent the evolution according to (3.16) of the 
initial operators Pj in (6.5), up to the time ti and ^2, respectively (see Appendix C). 

The two P-mesons that come from a decay of an T are quantum-mechanically entan- 
gled, in a way very similar to that of two spin 1/2 particles coming from a singlet state. 
As in that case, correlated measures on the two particles become physically significant. 
Indeed, the typical observables that can be studied at P-factories are double decay rates, 
i.e. the probabilities ^(/i, ti; /a, ^2) that a meson decays into a final state /i at proper 
time ti, while the other meson decays into the final state /a at proper time t2. They can 
be computed using: 



(6.7) 



g{fuh-j2,t2)^Tl [Of,® Of,) PAihM 



where C/^, Of, represent 2x2 hermitian matrices describing the decay of a single meson 
into the final states /i, /2, respectively. They can be identified with one of the two matrices 



^ Although other possibilities are in principle conceivable, this choice is the most natural 
one: it assures a consistent single meson dynamics when tracing over the degrees of freedom 
of one of the two P-mesons. [23] 



19 



in (4.6) and (4.8). Let us denote with O^, i = 1,2, 3, 4, the entries of the matrix Of when 
projected along the operators P^: 



The observables (6.7) can then be rewritten as 

g{h,h;h,h) = ^r-^ E n.(ti,t2) , (6.9) 



2 |r5 + tl 

where the elementary probabilities 



Py(ti,t2)=Tr{PiPi(ti)} Tr{P,P2(t2)} + Tr{P, ^2(^1)} Tr{P, ^(ts)} 

-TT{P,Ps{t^)}Tr{PjP^{t2)} - Tr{PiP^{t^)}Tr{PjPs{t2)} , 



(6.10) 



can be easily computed to any given accuracy using the solution of the evolution equation 
(3.16) collected in Appendix B. 

The expressions for the double decay rates in (6.9) can be compared with the results of 
the experiment. However, much of the analysis at P-factories is carried out using integrated 
distributions at fixed time interval t = ti — 12] this is a consequence of the short P-mesons 
lifetime and rapid B^-B^ oscillations, that does not allow a precise enough study of the 
double time dependence in (6.9). One is then forced to construct single-time distributions, 
defined by 



/•oo 

P(/i,/2;t)= / dt'g{h,t' + t;f2,t') 
Jo 

p. V o) n,Jt) , 



(6.11) 



12 



where 



/"OO 

Pi, (t) = / dt' V,jit' + t,t') , (6.12) 
Jo 

and t is taken to be positive. For negative t, one defines: 

/•oo 

r(/i,/2;-|t|)= / dt'g{h,t'-\t\;f2,t')e{t'-\t\); (6.13) 

the presence of the step-function is necessary since the evolution is of semigroup type, with 
forward in time propagation, starting from zero (we can not propagate a P-meson before 
it is created in a T-decay). In this case, one easily finds: 

r(/i,/2;-|t|) = P(/2,/i;|t|) . (6.14) 
20 



In the following, we shall always assume: t >0. 

The explicit form of the single-time probabilities nij{t) in (6.12) can be found in 
Appendix C; inserting these in (6.11), one is then able to compute the evolution in t of the 
double decay rate i^(/i, f2', t) for any specific decay state fi and /2. Before presenting the 
results for some experimentally relevant cases, let us note that in general the non-standard 
dynamics in (3.16) gives results for /^(/i, /2; t) that are quite different from those obtained 
in the usual case {i.e. in absence of JC). The most striking difference arises when the final 
states coincide /i = /2 = / and t approaches zero. Due to the antisymmetry of the initial 
state \iI^a) in (6.1), quantum mechanics predicts a vanishing value for i^(/, /; 0), while in 
general this is not the case for the completely positive dynamics generated by (3.1). This 
explains why correlated mesons systems turn out to be the most natural place to look for 
effects leading to loss of phase coherence and dissipation. 

The quantities i^(/, /;t), with t small, are therefore very sensitive to the dissipative 
parameters in (3.22) . This can be explicitly shown by considering the following combination 
of decay rates involving semileptonic final states: 

. rih+,h+;T) + rih-,h-:T) 

= rih^,h-,r)+r{h-,h^;r) ' ^'-''^ 

where again the variable t = tV has been introduced. Neglecting terms containing 5 times 
small parameters, one finds: 

7^(r) 





(1 + 


\rs 


2 




2) ( cosh 5t — cos wr) 


{\rs 


2 + 1 




2) cosh 6t + 2 TZe{rsr1) coswr 



+ 



(coswT + ly 



(6.16) 

I 

1 — iijo J 



As r approaches zero, the first term becomes vanishingly small, while the surviving pieces 
are all proportional to the dissipative parameters, so that: 

7e(0) = ^|A+^^[a;Xm(S)-MS)]| . (6.17) 

A non vanishing result in the measure of this quantity would clearly signal the presence 
of dissipation. However, to obtain more detailed information on the magnitude of the 
non-standard parameters, one has to study the full r-dependence in (6.16), and further 
analyze the behaviour of other observables that can be constructed out of the double 
probabilities F . 

For instance, one can reverse the sign in the numerator of (6.15) and study the asym- 
metry: 

_ r(fe+,fe+;r)-r(fe-,fe-;r) 
^^^^^ - r{h+,h-,r) + r{h-,h+;r) ' ^^'^^^ 

Unfortunately, the same phenomenon already noticed in the case of single meson asymme- 
tries occurs also here: the dependence on the dissipative parameters is always suppressed 
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by a factor 5, so that it becomes extremely difficult to extract information about the 
non-standard effects from these observables. 

The situation is different for the asymmetry ^Am(^)) analogous to the one in (5.10) 
for single-meson systems: 

[rih+,h-,T) + rih-,h+;T)] - [rih+,h+;T) + rih-,h-,T)] 
^Am[r) r{h+,h+;T) + r{h-,h-,T) + r{h+,h-,T) + r{h-,h+;T) " ^ ^ 

In fact, even in the case 6—0, this asymmetry has a very distinctive dependence on the 
non-standard parameters: 

AAmi'!') = 7 coswT-l- r7^e(6')l (l— COS lot) + (u> COS UT+ sin ut) Tie [ — ] ■ (6.20) 

\1 — iu) J 

Thanks to the different time behaviour, an experimental study of this observable should 
allow to extract precise information on A and the combination ujXm{B) — 7le{B); the 
analysis of TZ{t) in (6.16) should then provide an estimate on the parameter D. 

The asymmetry in (6.19) is presently under intensive study at S-factories: [51-53] it 
is used to obtain a precise determination of the mass difference Am, since in absence of 
dissipation: A Ami'!') — cos lot. Although an accurate analysis is certainly needed in order 
to estimate precisely the sensitivity of those experiments to the different r-dependent 
functions in (6.20), on the basis of estimated data acquisition and relative errors, one can 
reasonably expect that the actual measured data will at the end constrain the dissipative 
constants A and LoXm{B) — TZe{B) to a few percent level. This would allow to establish 
significant limits on the presence of dissipative effects in the B^-B^ system. 

To get information on the magnitude of the fourth parameter C in (5.6), one has 
to study observables that involve also hadronic decay channels, in particular those with 
definite CP-parity. A typical asymmetry that is measured at S-factories involves double 
decay rates P(/i,/2;t) in which one of the final states is semileptonic and the other is 
hadronic: 

r(;. + ./;r)-r(/,-./:r) 
■^'("^ = r(A+./;r) + r(/.-,/;x) ' '^''^'^ 

Since the two mesons are quantum mechanically correlated, neglecting AB = AQ violating 
effects, the semileptonic decay effectively "tags" the flavour of the meson that decays into 
/. In view of this, one usually consider the asymmetry in (6.21) and the one in (5.14) 
for single meson decays as the same observable; [1, 2] for instance, when / = J/ipKs both 
asymmetries can be directly expressible in terms of one of the angles parametrizing CP 
violation in the standard model. The situation radically changes in presence of dissipation; 
as already stressed, the form of the double decay rates r{fi, /2; r) turns out to be very dif- 
ferent from that predicted by ordinary quantum mechanics and this makes the asymmetry 
(6.21) distinct from that in (5.14). 

In presenting the explicit expression of A fir), one can set 5 — 0, since this choice 
does not affect the dependence on the dissipative parameters. Further, we shall assume 
that CPT-violation in "mixing" be small, and therefore neglect terms containing powers of 



22 



the parameter TZe{$) higher then one. Within these approximations, one exphcitly finds: 



^2Jm(A/) 2C/ 



+ 



+ 



1 + |A/|2 u "'\2 + iuj 
1-|A/|2 27^e(A/) 



sm OJT 



(6.22) 



1 + |A/|2 1 + lA/P 



ne{9) 



2C/ 
a; 



Tie 



2i + u 
2 + 100 



C 



cos UJT , 



where the decay parameter A/ is defined in (5.15), while Q represents again the intrinsic 
CP parity of the state /. 

Assuming the vahdity of the AB = AQ rule and neglecting CPT- violating effects in 
"mixing" = 0) as well as "direct" CP-violations (|A/| = 1), a fit of (6.22) with experi- 
mental data allows the determination of both TZe{C) and Xm{C), together with Xm{\f). 
The asymmetry ^/ , in particular for the final state f — J/ip Ks , will be measured with 
increasing accuracy at S-factories, so that a high sensitivity on the different r dependences 
in (6.22) is expected. [51] Therefore, the study oi Af could result in one of the best tests 
on the presence of dissipative effects in B physics. 

As remarked before, the form of «4/(t) in (6.22) is valid only for r > 0. For negative 
times, one has to use the an analogous expression, obtained from the definition in (6.21) 
by exchanging the positions of the semileptonic and / final states in the probabilities F. 
Explicit computations gives: 



M-V\) =^^#T^ ^e(^) + 7^e(A, - + 2 C/ y.) - ^ 7^ef C 



+ |A/|^ 



2 + iuj 



+ 



UJ 

27^e(A/) 



1 



smu; T 



(6.23) 



|Ari2 



ne{e) 



2C/ 



UJ 



Xm{C) 



cosa; T 



Note that the two expressions (6.22) and (6.23) coincide at r = 0, while in absence of 
dissipation (C = 0), (6.23) can be obtained from (6.22) by letting r — > — r; the fact that 
this no longer true for C 7^ is a clear sign of presence of irreversibility. 

The observables involving correlated S-mesons that have been discussed so far are 
accessible at the so-called asymmetric S-factories, where the T-decays take place in a 
boosted reference frame with respect to the laboratory. This allows the subsequent decays 
of the two neutral B mesons to be physically separated and therefore a good determination 
of the time-difference r. In symmetric -B-factories, this measure is impossible: the two 
mesons decay too quickly to allow a complete reconstruction of the events. One has 
then to resort to time-independent probabilities, obtained by further integrating the rates 
nfij2;t) in (6.11): 



n/1,/2) 



(6.24) 
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The asymmetries that can be constructed out of these quantities are not very sensitive 
to the non-standard parameters, because of the suppression by factors 5. Nevertheless, one 
can still form useful observables involving semileptonic final states; these are the total -B^- 
mixing probability: 

^ ^ nh+,h+)+r{h-,h-) 

r{h+,h+) + r{h-,h-) + r{h+,h-) + r{h-,h+) ' ^ ' ' 

and the ratio of the total, same-sign to opposite-sign semileptonic rates: 

rih+,h+)+r{h-,h-) 
- r{h+,h-) + r{h-,h+) ■ ^^-^^^ 

Both observables have a quadratic dependence on S; the approximation 5 = is therefore 
very accurate. With this choice, one explicitly finds: 



(1 - [7^e(^)]2) 4 



where the dependence on the non-standard, dissipative parameters occurs via the combi- 
nation: 

X = A + oj'^ D + uj[ujXm{B) -Tie{B)\ . (6.28) 

Independent measures of these two quantities would provide a way to estimate both 
\TZe{9y\^ and X, and therefore give limits on both dissipative and CPT-violating effects. 
Unfortunately, both observables are not very well known; Xb is the better determined 
parameter and the most recent data'*' give: xb = 0.198 ± 0.019. [54] Assuming [JZe{6)]'^ to 
be negligible and using the world average for a;, [55] from (6.27a) one gets the estimate: 
X = (5.7±4.8)xl0-2.t The accuracy on the determination of X will be greatly improved 
when also the measurements at the asymmetric S-factories will be available. Indeed, 
preliminary estimates of about one percent sensitivity to Rb have been reported to be 
attainable in [38], while a conservative estimate of two percent accuracy in the measure 
of Xb is indicated in [40]. If confirmed by the actual data, these sensitivities will allow 
a determination of the combination X with about a few percent accuracy, providing an 
interesting test on the presence of non-standard effects in S-physics. 



The determination of xb in [54] actually involves so-called semileptonic- type decays 
[36] rather than pure semileptonic final states; one can check that the form (6.27a) for the 
observable xb is valid also in that case. 

^ When a = 0, the inequalities (2.13) further imply a = 7 and b = c = (3 = 0, or 
equivalently A = B = D, C — 0. In this simplified case, an estimate on the variable X 
translates into a corresponding one for the surviving dissipative parameter a; in particular, 
from the above quoted value for X, one has: a = (4.5 ± 3.8) x lO"-*^^ GeV. 
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7. DISCUSSION 



The description of open quantum systems in terms of quantum dynamical semigroups 
allows a general and consistent approach to the study of physical phenomena leading 
to irreversibility and dissipation. When applied to the analysis of the propagation and 
decay of neutral meson systems, it gives precise predictions on the behaviour of relevant 
physical observables: the new, dissipative phenomena manifest themselves through a set 
of phenomenological parameters, a, 6, c, a, (3 and 7, whose presence can be experimentally 
probed. 

Indeed, as discussed at length in the previous sections, various observables involving 
B^-W decays can be identified as being particularly sensitive to the new, non-standard 
effects. These observables will be measured with great accuracy in the new generation 
of dedicated fi-experiments, both at colliders (CDF-II, HERA-B, BTeV, LRC-b) and at 
5-factories (BaBar, Belle, CLEO-III), so that stringent bounds on the dissipative effects 
can be expected in the future. 

From the experimental point of view, the actual visibility of these effects clearly de- 
pends on the magnitude of the parameters a, b, c, a, (5 and 7. A pure phenomenological 
approach can not provide such information. However, in the framework of open quantum 
systems, the neutral 5-mesons are described as subsystems in interaction with an envi- 
ronment. In such instances, the effects of irreversibility and dissipation can be roughly 
estimated to be proportional to the square of the S-meson mass divided by the character- 
istic energy scale of the environment. Assimilating this scale to the Planck mass results 
in a very small estimate for the magnitude of the dissipative parameters, roughly of order 
10~^^ GeV. However, the sophistication of the experiments mentioned before is so high, 
that the sensitivity needed to provide useful constraints on such tiny effects should be 
reached in just a few years of data taking. [49-53] 

Being based on the general theory of open systems, the above estimate on the magni- 
tude of the non-standard, dissipative effects is rather robust and quite independent from 
the details of the actual dynamics that drives the interaction between subsystem and envi- 
ronment. Nevertheless, it has been recently questioned, [56] on the basis of the similarity of 
a simplified version of the evolution equation in (3.1) with those describing the dynamical 
reduction of the wave-packet. We point out that this analogy is only superficial: the phys- 
ical effects leading to dissipation in open systems are clearly distinct from those advocated 
as responsible for the reduction process; and indeed, the quantum dynamical semigroup 
generated by (3.1) can not describe dynamical reduction phenomena. 

The arguments of [56] have been further used in [57] to support the claim that non- 
linear evolution equations should be used to analyze dissipative effects in neutral meson 
systems, although also there dynamics of the form (2.2) are nevertheless adopted at the 
end. In the framework of open systems, the situation can be easily clarified. 

The dynamics of a small system S in interaction with a large environment E is in 
general very complex and can not be described in terms of evolution equations that are 
local in time: possible initial correlations and the continuous exchange of energy as well as 
entropy between the S and E produce memory effects and non-linear phenomena. Never- 
theless, when the typical time scale in the evolution of the subsystem S is much larger than 
the characteristic time correlations in the environment, the subdynamics simplifies and a 
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mathematically precise description in terms of quantum dynamical semigroups naturally 
emerges. [4-6] 

This description is very general and is applicable to all physical situations for which 
the interaction between S and E can be considered to be weak and for times for which 
non-linear disturbances due to possible initial correlations have disappeared. [7] These are 
precisely the conditions that are expected to be fulfilled in neutral meson systems: the char- 
acteristic time correlations in the environment, induced by the fundamental {e.g. gravita- 
tional or "stringy") dynamics, is certainly much smaller than the neutral meson lifetime, 
and the interaction between mesons and environment is for sure weak (its effects have 
not yet been detected). Furthermore, this open system paradigm automatically assures 
the fulfillment of basic physical properties, as forward in time composition and entropy 
increase (irreversibility). Therefore, the physical motivations for adopting a quantum dy- 
namical semigroup description of the extended neutral-meson effective dynamics appear 
to be rather compelling and general. 

As a final remark, let us point out that once effective dynamics generated by equations 
of the form (2.2) are accepted, the condition of complete positivity is then absolutely 
necessary for a physically consistent description of the neutral meson system. Indeed, 
time evolutions that do not satisfy this property unavoidably gives unphysical results. [42] 
This can be easily shown by examining the sign of the eigenvalues of of the density matrix 
pU(t, t) representing a correlated meson state. In fact, one of the properties that any density 
matrix needs to satisfy is that its eigenvalues be non-negative, for all times; without this 
basic requirement, its standard probability interpretation would be meaningless. 

As an example, following [28, 19, 20], let us consider a dynamics for pA generated 
by single-meson evolution equation of the form (3.1), where C is as in (2.14), but with 
vanishing a, b and c. The conditions of complete positivity in (2.13) are then clearly 
violated (unless a = 7 and /3 = 0). For 7i = 0, the eigenvalues Xi{t), i — 1,2,3,4, of 
the 4x4 matrix pA{t,t) can be explicitly obtained; being physical quantities, they can be 
computed in any phase convention, and in particular in the one for which the dissipative 
parameters in (3.22) take the simplified form A — B — a, C — i(3 and D — 'j. Then, one 
has: 

Ai,2 (t) = 2 ± { [E+ (t)] ' + [E_ (t)] % 2 [F(t)] n , 

1/2 (7-1) 
XsAt) = ±[E+{t) + E_{t)] { [E+{t) - E_{t)]'+4[F{t)]'] , 



where 



E±{t) = + 2a) * - (iv_ + 2a) e''^ *] , 



and v± = — (a -|- 7) ± A/(c^-^^yp~T^^ are both negative due to the simple positivity 
condition q;7 > As clear from the previous formulas, for t 7^ two of the eigenvalues of 
PA{t,t) are always negative, signaling the presence of unphysical "negative probabilities". 
When complete positivity is enforced though, one has = v- and all eigenvalues in (7.1) 
remain non-negative for all times. 
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These results offer a further motivation for studying generahzed dynamics of the form 
(3.1) at meson factories. These set-ups arc in fact high-performance quantum interferom- 
eters: at least in principle, they can clarify from the experimental point of view the role 
of the condition of complete positivity in the time evolution of correlated mesons. [42] As 
shown by the simple example above, this is not just a mere technical question: it is crucial 
to our physical understanding of the quantum dynamics of open systems. 



APPENDIX A 

As explained in the text, in order to find explicit solutions for the evolution equation 
(3.1), it is convenient to make a change of basis as described in (3.14). The resulting dy- 
namical equation has a simplified hamiltonian part TYq , while the dissipative piece becomes 
more involved: C = V JCV~^ . Its entries can be expressed as linear combinations of the 
dissipative parameters a, b, c, a, j3 and 7. One explicitly finds: 



{A.l) 



\rs + rL\ 



- A 


E 


A 


A*- 




A 








-A* 


n 


e 


. 


-A 


Q* 





where 



A = a(|rs|' - 1) (|rLp - l) + 27^e{(c - ib) [{{rL^ - l)rs - {\rsf - l)rL\ } 

-2(a + 7)7^e(rsr£) +2 7^e|(a-7 + 2^/3)rsrL} , 

E = a{\rL\^ - if - 47^e|(c - ib){\rL\'^ - I^l] + 2{a + ^)\rL\^ 

-27^e{(a-7 + 2^/3)ri} , 

S = a{\rs\^ - if + 47^e|(c - ib){\rs\'^ - l)rs} + 2{a + ^)\rs\'^ 

-27ee{(a-7 + 2z/3)r|} , 



A = a(|rL|^ - 1) [rsrl + l) + (a + 7)(^s - rL)rl + (c - ib){rs - tl - 2rs\rL\'^) 
+ (c + ib) [{tl - rs)rl - 2] + (a - 7 - 2ip)rl^ - (a - 7 + 2if3)rsrL , 
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$ = a{\rsf - 1) {rsrl + l) + (a + -f){rl - rj)rs + (c - ib) [{rl - r*s)rs + 2]rs 

+ (c + ib) {r*s -rl + 2\rs\^rl) - (a - 7 - 2i/3)r^r£ + (a - 7 + 2i/3)r| , 



e = -a(r5rL + l) - 2(c - z6) (r^rL + l)rL + 2(c + ib) (r^rL + l)r5 + 2{a + 7)r5rL 

+ (a - 7 - 2z/?)rj2 + (a - 7 + 2z/?)ri , 



1] 



-a|rsr2 + l|' + 27^e{(c - [(|r5p - l)rL - (Ir^p - l)rs] } 

- (a + 7)(ks|^ + I^lH - 27ee|(a - 7 + 2i/3)rsrL} 



Recalling the discussion in Section 3, an independent phase change of the basis vectors, 
\B^) ^ e'"^ I^O), 1^0) e"^ induces a transformat ion on the entries of the matrix C 
in (2.14), realized by the operator U^p in (3.6); explicitly, one finds: 

(c - ib) e-^^'t'-^^ (c-ib) , (a - 7 - 2z/3) ^ e^*^'^-^) (a - 7 - 2z/3) , {A.2) 

while a and a + 7 remain unchanged. Similarly, also the two ratios rs and in (2.9) 
change, according to the rule: 



rs e*^"^-^) rs , tl ^ e'^'^'^^ tl 



(A.3) 



As a result, the entries of the matrix C listed above are manifestly rephasing invariant. 



APPENDIX B 

In order to study the time evolution of B^-B^ observables, one has to solve the evo- 
lution equation (3.16), 



(B.l) 



with C as in (3.21), for a given initial state |p(0)). In other terms, one has to compute 
the entries of the 4x4 evolution matrix Mij{t), which gives the components pi{t), P2{t), 
P3{t)i P4{t) of the state vector \p{t)) at time t, in terms of the initial ones at t = 0: 



p,{t) = J2Mij{t)pj{0) 



i = 1,2,3,4 . 



(5.2) 
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As explained in the text, it is sufficient to write down an approximate expression for Mij{t) 
that contains contributions up to first order in the dissipative parameters (3.22) appearing 
in C. The expansion of Mij{t) within this approximation can be conveniently organized 
as the sum of two contributions: 



The matrix M^^^ has only diagonal non- vanishing terms: 

The entries of M^^^ take instead the following explicit expression: 



(S.3) 



(S.4) 



2C* 
Afl 



2C 



2C* 

2^* 
Afl 

2Am 



(S.5) 



Mi;)(t) 



Ar_ ' 

iB* 
2Am 



Mil^(t) = 0, 



2C 
ATV 



where Ar± = AT ± 2zAm. In presenting the above expressions for the entries of M^-^') 
and M(i) we have reconstructed the exponential dependences out of first order correction 
terms which are linear in time. [21] 

Although the time dependence in the solution of (-B.1) is always exponential, this 
feature is lost in perturbation theory. However, to a given order in the perturbative 



29 



expansion, one can always reconstruct the correct exponential behaviour by a redefinition of 
the parameters 75, 7/, and Am, so that they coincide with the widths and mass difference of 
the physical states \Bs) and |-Bl). To first order, only the widths get shifted: 75 ■ys + D, 
"iL ^ IL -\- D. Accordingly, if we redefine V to be the average of the new 75 and 71,, then 
the quantities F-i- get changed as: T±^T± + A — D. This explains the form of the 
exponential terms in {BA). 



APPENDIX C 



As explained in the text, the evolution in time of the two correlated neutral S-mesons 
coming from the decay of the T(45') resonance can be obtained using the single meson 
dynamics |p(0)) — > \p{t)) generated by the equation (3.16). It results convenient to regroup 
the four components of the vector \p{t)) into a 2 x 2 matrix 



~P{t) 



Pi{t) Pz{t) 



it can be expressed in terms of the elementary operators Pj, i 
(6.5): 



~p{t) = Y,PiPiif) ■ 



(C.l) 

1,2,3,4, introduced in 
(C.2) 



In order to compute explicitly the form of the various observables involving correlated 
mesons, we need to determine the time evolution of the operators Pi. Using the results of 
Appendix B, on easily finds [adopting the same notation as in (C.l)]: 



M\V{t) 



13 it) 

(1) 
43 



M. 



Mg\t) 



P2{t) = 



P4{t) = 



(C.3) 



M- 



1)/ 
34 

(1)/ 



(C.4) 



With these expressions, on can now calculate the elementary double probabilities in (6.10) 

n,(ti,t2) = Tr{P,P^(ti)} Tr{P,P2(t2)} + Tr{P^ P2(ti)} Tr{P^- P^(t2)} 

-lV{PiP3(tl)} Tr{P,P4(t2)} - TT{Pih{ti)]TT{Pjh{t2)] , 
and, with a further time integration, the single-time probabilities 

/•oo 

i7y(0= / dt' Vij{t' + t,t') . (as) 

Introducing again the variable t = Tt and recalling the definitions (5.3) and (5.6), one 
explicitly finds: 
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i7n 

i7l2 

ill3 
i7i4 

i72i 
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il23 
i724 
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e-^ D 
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2T ' 



cosh 5t + 



sinh 5t 



2T 5 + iuj 



i7i3(T) 
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sinh St 



2T 5-iuj 
il23(r 
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2T 5-%u 
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e-^ 


2T 1 


il34 
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e~'^ e 


2T 
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[t) = 


"i73i(r) 


il42 


[t) = 


i732(T) 


i743 


[t) = 


i734(r) 


il44 


[t) = 


i733(T) 



To; cos u>T + sin u>t] , 



where the parameter r represents the following ratio: 

2-5 + ioj 



(C.6) 



(C.7) 



2 + 6-iu 

These expressions for the components of nij{r) have been used to compute the double 
decay rates -r(/i,/2;T) in (6.11). 
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